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Abstract 


IS-95 is the only CDMA based second generation standard for cellular communication. 
Its performance has been evaluated over Rayleigh fading channel by many investigar- 
tors. However, Nakagami-m channel model is considered to be more appropriate for 
frequency selective fading channels. It includes the Rayleigh fading channel and one 
sided Gaussian as special cases. In this thesis the BER performance of forward link 
(base station to mobile station) of IS-95 standard has been evaluated. The multipath 
channel model assumes independent paths with Nakagami-m fading statistics. RAKE 
receiver has been used to evaluate the performance under various multipath fading 
conditions. The results indicate that for m<l, the performance is worse when com- 
pared to the performance with Rayleigh fading channel model, whereas it gets better 
for m>l. Also, more than 20 users can access the channel (assuming 64 chips Walsh 
code) simultaneously and can achieve uncoded BER of 0.01 with an average SNR of 
around 10 dB. 



Contents 


1 Introduction 1 

1.1 Spread Spectrum Systems 1 

1.1.1 DS-CDMA 2 

1.2 Channel Considerations 4 

1.3 Motivation for the Present Work 4 

1.4 Organization of the Thesis 5 

2 Review of Fading Channels 6 

2.1 Features of fading channels 6 

2.1.1 Slow and Fast fading 7 

2.1.2 Frequency-Flat and Frequency-Selective fading 7 

2.2 Modeling of fading channels 8 

2.2.1 Nakagami-m model 11 

2.2.2 Characteristics of Nakagami-m distribution 12 

2.3 Diversity Techniques 13 

2.3.1 Diversity combining techniques 14 

3 Downlink Analysis ^7 

3.1 Probability of Error for Nakagami-m channel 20 

4 Simulation Methodology 22 

4.1 Transmitter 22 

4.2 Channel 25 

4.2.1 simulation steps for m > 1 25 

i 



4.2.2 simulation steps for m < 1 27 

4.3 Receiver 29 

5 Results and Conclusion 31 

5.1 Scope for Future work 36 

References 37 



List of Figures 


1.1 Block diagram of DS-CDM A transmitter 3 

1.2 Block diagram of DS-CDM A receiver 3 

2.1 Multipath propagation model[3] 9 

2.2 The distribution of delays[ll] 10 

2.3 Nakagami PDF for fi! = 1 and various values of the fading parameter m. 11 

2.4 Block diagram for EGC and selection combining techniques [11] .... 15 

2.5 Block diagram for MRC technique [11] 15 

3.1 Multiuser modulator and multipath. chan.nel[3] , , , , , •. •. , •. •. •. 15 

3.2 Rake demodulator for multipath propagation [3] 19 

4.1 Block diagram of transmitter for single-user [11] 23 

4.2 SSRG configuration for generation of (a) In-phase PN sequence (b) 

Quadrature PN sequence [11] 24 

4.3 FIR Filter specifications [11] 24 

4.4 Impulse response of FIR filter [11] 26 

4.5 Simulation model for Nakagami-m fading channel 28 

5.1 Autocorrelation function of I-phase PN code 32 

5.2 Autocorrelation function PN code of length 15 32 

5.3 Simulation results for m= 0.5, 0.7 and 1.0 33 

5.4 Simulation results for m= 1, 3, 5 and 9 34 

5.5 Results for m=3 with varying number of users 34 

5.6 Curves for Chernoff bound and simulation result for m=l 35 

5.7 Curves for Upper bound, Chernoff bound and simulation result for m=3 35 

111 



List of Tables 


4.1 CoeflScients of FIR filter [11] 


26 


IV 



Chapter 1 


Introduction 


In recent years, there has been considerable interest in the application of Code Divi- 
sion Multiple Access (CDMA) techniques to multiple access communications. CDMA 
can potentially accommodate more users than either Time Division Multiple Access 
(TDMA) or Frequency Division Multiple Access (FDMA), because as additional users 
are activated, the system performance as seen by any single user is degraded only s- 
lightly. Thus, as long as the increase in Bit Error Rate (BER) is tolerable, the already 
established links remain operational. Another advantage of CDMA is that there are no 
requirements for precise time or frequency coordination among users as in TDMA or 
FDMA. Furthermore, CDMA offers the usual benefits associated with spread spectrum 
systems such as interference rejection/suppression and multipath mitigation. 


1.1 Spread Spectrum Systems 

The spread spectrum system is a wideband system which expands the system band- 
width several times the message bandwidth. It has been defined as [1] Spread spec- 
trum is a means of transmission in which the signal occupies a bandwidth in excess of 
the minimum necessary to send the information; the band spread is accomplished by 
means of a code which is independent of the data, and a synchronized reception with 
the code at the receiver is used for De-spreading and subsequent data recovery.” The 
applications of spread spectrum technique are many fold[2, 1]: 
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1. Interference suppression 


2. Energy density reduction 

3. Ranging or time delay measurement 

4. Anti-jamming 

5. Multiple user random access communications with selective addressing capability 

6. Anti-interference 

7. Low probability of intercept 

8. Accurate universal timing 

A significant feature of a spread spectrum system is its Immunity to multipath recep- 
tion. In a wireless system a signal transmitted may reach the receiver via different 
paths ( being reflected and diffracted by various terrain configurations) apart from a 
direct path. Also there are situations when no direct path is available. This results in 
a distorted received signal, which is also known as self jamming. But a spread spec- 
trum system is resistant to such situations and diversity techniques can be employed 
to separate different component path signals which can be processed further to demod- 
ulate received signal. Because of the possibility of random access communication, the 
spread spectrum communication system is also called spread spectrum multiple access 
(SSMA) communication. The SSMA system is divided into Direct sequence CDMA 
(DS-CDMA), Erequency Hopping CDMA (FH-CDMA) and hybrid CDMA. DS-CDMA 
is most important amongst all three and hence it is only being considered in this thesis. 

1.1.1 DS-CDMA 

DS spread spectrum signals are generated by linearly modulating the data stream with 
a wideband signal. This wideband signal is known as signature sequence. This sequence 
is assigned to a user and is unique to each user. The block diagram of a DS-CDMA 
system is shown in Fig. 1.1. In a mobile communication environment, the important 
characteristics of the DS-CDMA are universal frequency reuse, eflScient power control. 
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soft hand off of the mobile, and enhancement of capacity using sectored antenna etc. 
[3]. Referring to Fig. 1.1, the data source is modulated by the PN sequence of the 



Transmitter 

Output 


Figure 1.1: Block diagram of DS-CDMA transmitter 



Decision 

Device 



Demodulated 

Output 


Figure 1.2: Block diagram of DS-CDMA receiver 

user followed by carrier signal. The PN sequence is a stream of I’s and O’s whose 
rate is manifold higher than the data rate. Individual bit of a PN sequence is called a 
chip. The ratio of the PN sequence rate to the data rate is known as processing gain 
of the system. The receiver structure for the DS-CDMA system is shown in Fig. 1.2. 
The received signal from the channel output is carrier demodulated and then passed 
through a correlator which correlates the carrier demodulated received signal with the 
locally generated PN code of the user. Here it is assumed that the locally generated 
PN sequence generator is in synchronism with that of the transmitter. The correlator 
output is fed to a decision device which makes some processing and detects the bits. 
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1.2 Channel Considerations 

A prerequisite for the performance analysis of mobile communication system is a model 
for channel behavior i.e. fading statistics. With an accurate channel model, it becomes 
feasible to compare various systems. Typical wave phenomena like scattering, diffrac- 
tion, and absorption by physical structures in the environment leads to a diffusion of 
a transmitted wave into a continuum of partial waves with different amplitudes and 
phases [7]. As a result, a transmitted signal, modulating a carrier wave is attenuated 
and distorted. These effects depend on the geometrical and electromagnetic properties 
of the physical transmission environment. Real environments show a broad variation of 
the physical properties, leading to varying channel behavior. Typically, the scattering 
structures are irregularly distributed. The scattered partial waves interfere spatially 
and build up an irregularly distributed interference pattern. Such circumstances have 
led to a large number of extensive experiments and numerous theoretical investigations 
to be performed on the intensity distribution of fading under various conditions. As a 
result diverse forms of the distributions have been presented. Amongst them, following 
three may be regarded, in view of practical uses, as the representatives [6]. 

• Rayleigh Distribution 

The applicability of this distribution to fading in many practical channels has 
been confirmed by many investigators. 

• Log-normal Distribution 

Log-normal law is often followed over long periods over which the short-term 
fading is nonstationary. 

• Nakagami-m Distribution 

This distribution degenerates to Rayleigh distribution for m = 1. 

1.3 Motivation for the Present Work 

The BER performance of each user in a CDMA system is affected by fading statistics 
of the received signal. A series of propagation experiments conducted in typical ur- 



CHAPTER 1. INTRODUCTION 


5 


ban/suburban areas have revealed the multipath nature of the mobile radio channel 
resulting from reflection, refraction, and scattering by buildings and other obstructions 
in the vicinity of the mobile [5, 14]. Multipath propagation can be resolved at the 
receiver of a DS-CDMA system since it uses a bandwidth much larger than the co- 
herence bandwidth of the channel. The problem of obtaining the distribution of the 
signal strength of each path coincides with the random phaser interference problem. 
Nakagami-m distribution provides a more general solution for such a problem as com- 
pared to Rayleigh and Rice distributions, which are only special case solutions [6, 7]. 
It has been shown that such multipath scenario in urban /suburban areas are better 
modeled by Nakagami-m distribution [4, 5, 6]. As a brief digression, the Nakagami-m 
distribution is also of great interest to ionospheric physicists because amplitude fading 
due to ionospheric scintillation has been found to follow an m-distribution [9]. 

It has been established that Nakagami-m distribution model is best fit for mobile 
communication system channel model. Though lot of work has been carried out to 
evaluate the performance of DS-CDMA systems over Rayleigh fading channel, the use 
of Nakagami-m channel model would be more appropriate. It would be of interest 
to analyse how IS-95, which is leading 2"*^ generation standard for mobile communi- 
cation, performs when incorporated through Nakagami-m channel. In this work, a 
channel model incorporating Nakagami multipath fading is developed and the BER 
performance of IS-95 standard using coherent demodulation has been evaluated. For 
simulation purpose only the forward link of IS-95 standard [11] has been used incor- 
porating the RAKE receiver [3] using Maximal Ratio Combining (MRC) technique 
[ 8 ], 


1.4 Organization of the Thesis 

Chapter 2 contains the details of the multipath channel characteristics and various 
methods to alleviate the effect of the interference resulting due to multipath. Chap- 
ter 3 gives detailed analysis of forward link of DS-CDMA system. In chapter 4 the 
implementation steps have been described in detail. Chapter 5 presents the results, 
conclusion and the scope for future work. 



Chapter 2 


Review of Fading Channels 


Any communication system has transmitter, receiver and in between the two is the 
channel. For mobile communication system the transmitter and the receiver may be 
mobile. When the data or voice is transmitted from the base station (BS) the trans- 
mitter is fixed whereas the receiver ( mobile station MS) is movable. The link from the 
BS to MS is called forward link and the link from MS to BS is termed as reverse link. 
Forward link communication is one to many i.e. there is one transmitter and many 
receivers and the reverse link communication is many to one i.e. one receiver receives 
signals from many transmitters. The channel of a mobile communication system has 
a tendency to change with time i.e. the impulse response of these channels are time- 
variant, this is the consequence of the constantly changing physical characteristics of 
the media. Also the transmitted signal is reflected and refracted by a variety of smooth 
and rough surfaces, so that it is replicated at the receiver with several time delays. This 
is called multipath propagation. 


2.1 Features of fading channels 

If we transmit an extremely short pulse over a time-varying multipath channel, the 
received signal will appear as a train of pulses. Hence one characteristic of a multipath 
medium is the time spread introduced in the signal that is transmitted through the 
channel. Now let us consider the transmission of an unmodulated carrier at frequency 
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/c) then the low pass equivalent of the received signal is given by[12] 

ri{t) = (2-1) 

n n 

Thus, the received signal consists of the sum of a number of time-variant vectors 
having amplitudes Q:„(t) and phases The fading phenomenon is a result of the 

time variation in the phase On{t) i-e- the randomly time varying phase at times results 
in the vectors adding destructively. The result is that the received signal is very small 
or practically zero. Thus, the amplitude variation in the received signal is termed as 
signal fading and is due to the multipath characteristic of the channel. 

2.1.1 Slow and Fast fading 

The notion of slow and fast fading is related to the coherence time Tc of the channel, 
which measures the period of time over which the fading process is correlated. The 
coherence time Tc can be defined as [13], the period of time after which the correlation 
function of two samples of the channel response taken at the same frequency but 
different time instants drops below a certain predetermined value. The coherence time 
is also related to the channel Doppler spread fi by 

Tc ~ ^ (2.2) 

The fading is said to be slow if the symbol time duration T, is smaller than the channel’s 
coherence time Tcj otherwise it is considered to be fast. In slow fading a particular 
fade level will affect many successive symbols, which leads to burst errors, whereas in 
fast fading the fading decorrelates from symbol to symbol. 

2.1.2 Frequency-Flat and Frequency-Selective fading 

If all the spectral components of the transmitted signals are affected in a similar man- 
ner, then the fading is said to be frequency nonselective or frequency flat. This is the 
case for narrowband systems in which the transmitted signal bandwidth is much small- 
er than the channel coherence bandwidth Wc- This bandwidth measures the frequency 
range over which the fading process is correlated. It is defined [13] as the frequency 
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bandwidth over which the correlation function of two samples of the channel response 
taken at the same time but at different frequencies falls below a suitable value. The 
coherence bandwidth is related to the TnoMfrium delay spread Tmax by 

Wa (2.3) 

‘^max 

Whereas, if the spectral components of the transmitted signals are affected by different 
amplitude gains and phase shifts, the fading is said to be frequency selective. This 
applies to wideband systems in which the transmitted bandwidth is bigger than the 
channel’s coherence bandwidth. 


2.2 Modeling of fading channels 


When fading affects mobile communication systems, the received carrier amplitude is 
modulated by the fading amplitude a, where a is a random variable (RV) with mean- 
square value and probability density function (PDF) Pa(o:), which is dependent 

on the nature of the radio propagation environment. After passing through the fading 
channel, the signal is affected by additive white Gaussian noise (AWGN), which is 
assumed to be independent of the fading amplitude a and is characterized by one- 
sided power spectral density No (W/Hz). Hence the received signal power is modulated 
by a^. Thus, the instantaneous signal-to-noise power ratio (SNR) per symbol can be 
defined as 7 = a^Ej,/No and the average SNR per symbol by 7 = QE^/Ng, where Ea 
is the energy per symbol. 

The PDF of 7 is given by [13] 


P'ri'r) = 


Pa{\/^lh) 

2^77/^ 


(2.4) 


The moment generating function (MGF) M^{s) associated with the fading PDF Py{y) 
is defined as 


rOQ 

Mj{s)- / p^{'y)e^'^d'y 

J 0 


(2.5) 
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The amount of fading (AF) or fading figure associated with the fading PDF is defined 
as 


^ p _ varja'^) _ £;({a* - fi)^] ^ E(j^} - 

When wideband signals propagate through a frequency-selective channel, their spec- 
trum is aflfected by the channel transfer function, resulting in a time dispersion of the 
waveform. This type of fading can be modeled as a linear filter characterized by the 
following complex valued lowpass equivalent impulse response 





x^(t> 


L Paths with 
varying ampiitude, 
phase and delay 


Figure 2.1: Multipath propagation model[3] 


f'P 

h(t) exp(-j6i)5(t - n) (2.7) 

i=i 

Where (J(.) is the Dirac delta function, I the path index, and . and 

are the random channel amplitudes, phases, and delays, respectively. The model for 
multipath is shown in Fig. 2.1. 

In Eqn. (2.7) Lp is the number of resolvable paths (the first path being the reference 
path whose delay =0). Obviously Lp is related to the ratio of maximum delay spread 
to the symbol time. Under the slow-fading assumption, Lp is assumed to be constant 
over a certain period of time, and and are all constant over a 
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Figure 2.2: The distribution of delays[ll] 

symbol interval. The fading amplitude oii of the resolved path is assumed to be a 
RV whose PDF be Nakagami-m distributed. 

With out loss of generality it can be assumed that the power component decreases 
with respect to delay, the first arriving path exhibits a lower amount of fading whereas 
the last arriving path exhibits highest amount of fading. The are related to the 

channel’s power delay profile (PDP), which is also referred to as the multipath intensity 
profile (MIP) and is decreasing function of delay. The mobile radio channel is well 
characterized by an exponentially decaying PDP 

Cll — 6Xp( 'ri/Tjnax)'i f ~ f > 2, , Lp (2-8) 

where Oi is the average fading power corresponding to the first (reference) propagation 
path and Tmax is the channel maximum delay spread. 

The delay spread is defined as the square root of the difference between the mean of 
the delay squared and the square of the mean delay, which is expressed as[ll] 

= v^f2 - f2 (2.9) 

The delays are often assumed to exponentially distributed [11] shifted by n as illus- 
trated by Fig. 2.2. If this distribution is normalized by its delay spread, then it is 
given by the expression 



r-Tj 

Ar 


) 


( 2 . 10 ) 
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The mean of the normalized exponential distribution is 

1 TOO r — Ti 

^)dr = A^ + ri (2.11) 

The mean square of the distribution is given by 

^ ^ exp( ■^)dT = A2 + (A^ + ri)2 (2.12) 

Thus the standard deviation for a delay with an exponential distribution is 

CTr = _ f2 ^ (2.13) 


2.2.1 Nakagami-m model 


The Nakagami-m PDF is in essence a central chi-square distribution given by [6] 


p(a) = 




exp(- 


), a > 0 


(2.14) 


where m is the Nakagami-m fading parameter which ranges from | to oo. Fig. 2.3 
shows the Nakagami-m PDF for = 1 and various values of m-parameter. 



Figure 2.3: Nakagami PDF for 0 = 1 and various values of the fading parameter m. 

Applying Eqn. (2.4) the SNR per symbol, 7, is distributed according to a gamma 
distribution given by 


^7(7) = 


1 

7”*r(m) 


exp( 



7 > 0 


(2.15) 
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The MGF in this case is given by 

M^(s) = (1 - (2.16) 

m 

and the moments are given by [ 6 ] 

which using Eqn. (2.6) yields an AF of 

AF = — , m > ^ (2.18) 

Hence, the Nakagami-rn distribution covers the widest range of fading figure (from 
0 to 2) among all the multipath distributions. It includes the one-sided Gaussian 
distribution (m = 5 ) and the Rayleigh distribution (m = 1 ) as special cases. In the 
limit as m — >■ + 00 , the Nakagami-m fading channel converges to a nonfading AWGN 
channel [13]. Finally, the Nakagami-m distribution often gives the best fit to land- 
mobile [5, 7] radio link. 

2.2.2 Characteristics of Nakagami-m distribution 

• If the envelope is Nakagami distributed, the corresponding instantaneous power 
is 

gamma distributed. 

• The parameter m is also called the ‘shape factor’ of the Nakagami distribution. 

• In the special case m = 1, Rayleigh fading is recovered, with an ‘exponentially 
distributed’ instantaneous power. 

• For m > 1 , the variance of random variable reduces as compared to Rayleigh 
fading. 

• Nakagami fading occurs for instance for multipath scattering with relatively large 
delay-time spreads, with different clusters of reflected waves. Within any one 
cluster, the phases of individual reflected waves are random, but the delay times 
are equal for all waves. As a result the envelope of each cumulated cluster signal 
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is Rayleigh distributed. The average time delay is assumed to differ significantly 
between the clusters. If the delay times also significantly exceed the bit time of a 
digital link, the different clusters produce severe intersymbol interference, so the 
multipath self-interference then approximates the case of co-channel interference 
by multiple incoherent Rayleigh-fading signals. 

• The Nakagami model is also often used to describe the interference cumulated 
from the multiple independently Rayleigh-fading sources, particularly if these are 
identically distributed. 

2.3 Diversity Techniques 

Errors in reception tend to occur when the channel is in deep fade i.e. when the channel 
attenuation is large. If we can supply to the receiver multiple replicas of the same 
information signal transmitted over independently fading channels, the probability that 
all the signal components will fade simultaneously is reduced. There are several ways 
in which we can provide the receiver with multiple independently fading replicas of 
the same information bearing signal. Depending on the propagation mechanism, these 
may include: 

• Space diversity 

We may employ a single transmitting antenna and multiple receiving antennas. 
Usually separation of at least 10 wavelengths is required between two antennas 
in order to obtain signals that fade independently [12]. 

• Frequency diversity 

The same information bearing signal is transmitted on different carriers, where 
the separation between successive carriers equals or exceeds the coherence band- 
width fc of the channel to obtain independently fading signals. 

• Angle of arrival diversity 

If more than one receiving beam oriented towards different non overlapping re- 
gions are employed, the signals received via these beams would be decorrelated. 
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• Polarization diversity 

rhis iriiplios a singlo polarization at the; transmitter, with depolarization in the 
propagation medium. Independent reception is possible with two orthogonal 
polarizations, and the two resulting signals do not fade in correlated manner. 

• Time diversity 

The same information bearing signal is transmitted in different time slots, where 
the separation between successive time slots equals or exceeds the coherence time 
Tc of the channel to ensure the independence of different fading signals. 

• Multipath diversity 

It can be obtained by using the signal having a bandwidth much greater than 
the coherence bandwidth of the channel. Such signal will resolve the multipath 
components and, thus, provide the receiver with several independently fading 
signal paths. If the signal bandwidth is W and multipath spread is of 7^ sec, 
there are TmW resolvable components. Since j^, fc being channel coherence 
bandwidth, the number of resolvable components may also be expressed as W/fc- 
Hence the order of diversity obtained by using this technique is Wj fc [12]. 

Of these, the last two have been primarily considered for digital data transmission. 
Except for polarization diversity, any number of more-or-less independently fluctuating 
signals can be made available by multiplication of equipment. 

2.3.1 Diversity combining techniques 

The common linear combining methods that are used are: 

• Selection combining 

This technique is based on choosing the channel with highest SNR and rejecting 
all others. Fig. 2.4 shows the schematic diagram for this technique. 

• Equal gain combining (EGG) 

In this case the decision statistics from all the L channels are combined with 
equal gain for final detection. The combiner is a post-detection combiner. Fig. 
2.4 gives the details for this technique. 
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X = E > for EGC 

1-1 


X — select the largest of L channels > for Selection diversity 


User 

Data 


Figure 2.4: Block diagram for EGC and selection combining techniques [11] 



User 

Data 


Figure 2.5: Block diagram for MRC technique [11] 
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• Maximal ratio combining (MRC) 

The phases of the diversity signals from L channels being added together are 
aligned and there envelopes are weighted in proportion to the square roots of their 
SNR’s [11]- This is the optimal combining technique as it yields the maximum 
achievable SNR. Block diagram for this technique is shown in Fig. 2.5 . 



Chapter 3 


Downlink Analysis 


When the spread spectrum employs'pseudo random sequence with chip time Tc, inverse- 
ly propoitional to the spreading bandwidth, the individual paths can be distinguished 
if they are mutually separated by delays greater than T^. Since in this case the various 
delayed versions of the signal will be mutually nearly uncorrelated. Received signal of 
A*'* user[3j 

____ _ r, 

JCi(t) — y £<c(^) ^n(^) ^ 5-^ Oi(h(t — nTc — 6i) X |a^^^(A:)cos [27r/o(i — 6[) + (f>j \ -|- 

n fel 

o!:^\k)sin [27r/o(t - (Ji) + 

where, Ec is the energy per chip of the data, Xn{k) is the nth symbol bit of the kth 
user, Qi is the attenuation of the 1th path, h{.) is the baseband wave shaping filter 
impulse response, 6i is the delay of the 1th path, (f>i is the phase of the 1th path signal, 
a^n^{k) and is the in phase and quadrature PN sequence of the user. 

Each of the distinguishable multipath components will actually be a linear com- 
bination of several indistinguishable multipaths of varying amplitudes. Since these 
will add as random vectors, the amplitude of each term will appear to be Nakagami-m 
distributed and the phase uniformly distributed. 

Let Ao(l) be the unmodulated pilot signal so that x„(0) = 1 for all n. We assume 
that the pilot pseudo random sequence is shared by all users. Hence the received signal 
containing users and a pilot sequence, all originating from the same Base station 
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Figure 3.1: Multiuser modulator and multipath channel[3] 
will be (Fig, 3.1) 

= + (3.2) 

k=:l 

Xo(t) is further scaled by Aq, the additional gain allotted the pilot signal. a^(0) and 
a^(0) are pilot’s QPSK spreading sequences. The users sequences are the products of 
those of the pilot and of the user-specific sequences b|^(k), b^{k) i.e. 

aUk) = aU0).bUk). «?(*) = «?(0).6?(A) (3-3) 

The timings of all individual users is locked to that of the pilot sequence, so that 
multipath delays need only be searched on the pilot sequence. 

The delays of different multipath components are estimated. Then it is used to 
remove pilot QPSK spreading, giving rise to the quadrature outputs 

yfK.[Ao + Xn{k).bl{k)].aiCOs{<l>i) + 

\[Ec\Aq + a;„(*).6^(fc)].oj5m(0i) + 


(3.4) 

(3.5) 
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From 

MatcNd 



Figure 3.2: Rake demodulator for multipath propagation [3] 


Ao is the pilot gain, and are the contribution of all other (uncorrelated) multipath 
components as well as those of all other users. From these outputs the relative path 
values aicos^i and aisin4>i can be estimated by simply averaging over an arbitrary 
number of chips Np. This number should be as large as possible without exceeding the 
period over which ai and (f>i remain relatively constant. 

The result of n** db ip of the path, after multiplication by the quadrature-user-specific 
pseudorandom sequences as shown in Fig. 3.2, is 

yin{k) = i/BcXn{k)QiaiC03{<^i--$i)+di{u^^cos^i + u^hin^i) (3.6) 

Summing over the N chips over which a:„(A:)is constant (±1), we have 

yi(ib) = ±N/EmdiC08i^i - it) + dilcosii ^ + sinii t'P] (3-7) 

Thus, 


= -11 = -Ar[v/Saiai<»»(A - A)l. 


(3.8) 
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with a change of sign if Xn{k) = l,and 


V'ar[Kt(*)] = Neil‘S 1q/2 


(3.9) 


where Iq is given by Iq = Nq + 1)*“ j Ec(j). 


3.1 Probability of Error for Nakagami-m channel 


Assuming that eac:h path’s attenuation is uncorrelated to that of others, the multipath 
component ainplituties be RV’s mutually independent, then the error probability for 
perfect estiinatP.s is given by [3] 

L 

n 

(3.10) 


,o,,])<Binexp(-afjVB,/4)J 


llt?re, E, NEc is the ^-cliip symbol energy, and expectations are with respect to the 
random variables ai. 

Taking ai variable to be Nakagami-m distributed the PDF of ai is given by (2.14), 
which i.s being rt'prtxliiced here 


PiOll) = exp(-— ), Oii > 0. 


r(m)nr fit 

where Qi = and letting A - «?. we obtain 


(3.11) 


pWi) 


_ 

dat 


-ypi 


exp(-mAM) g 


ftl"r(m) 

Thus, for Nakagami-m distributed attenuations, 

. E[e^] 


( 3 . 12 ) 


1 


Ti + S&F* 


(3.13) 
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Hence, the probability of error for L multipaths assuming all paths having equal average 
strength (i.e = fl, V 1) is 


(3.14) 


Where Eg — QEg. Eqn. (3.14) is the Chernoff bound on probability of error for 
Nakagami-m channel. 

Using MGF approach [13], the average BER for QPSK is given by 


where is given by (2.16) for Nakagami-m distribution, and ji — Hence, the 

average probability of error is given by the expression 




+ 


D.E. 


2mIoSi'n?(j> 


)-^^d4> 


(3.16) 


An upper bound on (3.16) can be obtained by evaluating the integrand a.t (p = tc/2, as 
integrand has a single maximum at this value. 


Fe< j(l + 


2mlg 


(3.17) 



Chapter 4 


Simulation Methodology 


This simulation has been carried out to evaluate the performance of forward link of 
DS-CDMA (IS-95 standard) over Nakagami-m channel. The assumptions made in the 
simulation are: 

• The delays of multipath do not vary significantly during the life of one call. 

• The fading parameters i.e. attenuation and phase remain nearly constant over 
one symbol duration. 

• The value of fading parameter m for Nakagami-m channel is assumed to be same 
for all the multipaths. 

• It has been assumed that the local oscillator at the receiver is synchronized with 
the transmitter. 

• The delay of each of the multipaths is known to receiver i.e. it has been assumed 
that the chip time of PN sequence has been accquired. 

4.1 Transmitter 

The transmitter of Fig. 3.1 has been implemented. The block diagram for IS-95 single- 
user transmitter is shown in Fig. 4.1. The implementation steps are as follows: 
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X Ch«nn«l PK cod* 

Chi.pii mt: 1.2288 Hcpm 

aamZni^ t 



Q Ch«zm«l PK codU* «i.n2 7ti^ t: 

Ch±p« at: 1.2288 Mcp« 


Figure 4.1: Block diagram of transmitter for single-user [11] 


• The input bit stream consisting of plus one’s and minus one’s (+1 and -1) has 
been generated using random number generator with uniform distribution. 


The Walsh function for 64 (N=32) chip length has been generated by [11] 




Hn -Hn 


The basic building block representing logic elements {0,1} as {1,-1} is 


1 1 
1 -1 


• The characteristic polynomials of in-phase and quadrature PN sequences are 
given by [11] 

//(a:) = 1+ -H r® + a:’’ + a:® + 

fQ{x) = 1 -f- a:® + a:^ -I- x® -f- a:® + a:^° + a;^^ + a:^^ + a;^® (4.3) 

These in-phase and quadrature sequences are generated by mechanizing the bi- 
nary primitive polynomials Eqn. (4.3) as simple shift register generators (SSRG) 
configuration as shown in Fig. 4.2. The initial vector of the registers in the figure 
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Clock 


(a) 



Clock 

(b) 


Figure 4.2: SSRG configuration for generation of (a) In-phase PN sequence (b) Quadra-- 
ture PN sequence [11] 



Figure 4.3: FIR Filter specifications [11] 
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is 000000000000001. The length of each PN sequence is 2^® — 1 = 32767. One 
zero bit is added at the end of the largest runlength of (14) zeros in order to make 
the length of PN sequences equal to 32768. In IS-95 this has been done to have 
exactly 75 cycles of the sequences in 2 secs, as the clock rate is 1.2288 Mcps. 

• The data bits axe modulated by the channel assigned Walsh sequence resulting 
64 chips (at the rate of 1.2288 Mcps) per bit. Then the Walsh coded data bits are 
modulated by I and Q phase PN sequence, with particular ojSfeet which identifies 
the base station, we get chips at the rate of 1.2288 Mcps at the output of PN 
sequence modulator. These chips are convolved with the coefficients of FIR filter 
to shape the baseband pulse. The filters are specified to have a frequency response 
H(J) that satisfies the limits shown in Fig. 4.3. The filter coefficients are are given 
in Table 4.1 [11] and impulse response is shown in Fig. 4.4. 

• The baseband waveform is then carrier modulated for QPSK modulation. 

4.2 Channel 

Simulation of the multipath channel model is carried out in accordaxice with the Fig. 

4.5. It has been assumed that each path constitutes several discrete paths, as a result 

the envelope of each path is Nakagami-m distributed [9]. 

4.2.1 simulation steps for m > 1 

• Select an appropriate value of the m parameter. 

• Generate m exponentially distributed random variables {Yi, Y 2 , , Ym), by first 

generating m uniform random variables (Xi,X 2 , ......... X^) using the random 

function, and then taking Yi = —log{Xi) for i = 1, 2, , m. 

• Average F = (li + F 2 + + Ym)/m, then Y is gamma distributed. 

• Take the square root ofY,Z = y/Y, where Z will be Nalcagami-m distributed. 



Amplitude 
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n 

ha(n) 

n 

/io(n) 

0, 47 

-0.025288315 

12, 35 

0.007874526 

1, 46 

-0.034167931 

13, 34 

0.84368728 

2, 45 

-0.035752323 

14, 33 

0.126869306 

3, 44 

-0.016733702 

15, 32 

0.094528345 

4, 43 

0.021602514 

16, 31 

-0.012839661 

5, 42 

0.064938487 

17, 30 

-0.143477208 | 

6, 41 

0.091002137 

18, 29 

-0.211829088 

7, 40 

0.081894974 

19, 28 

-0.140513128 

8, 39 

0.037071157 

20, 27 

0.094601918 

9, 38 

-0.021998074 

21, 26 

0.441387140 

10, 37 

-0.060716277 

22, 25 

0.785875640 

11, 36 

-0.051178658 

23, 24 

1.000000000 


Table 4.1: Coefficients of FIR filter [11] 



Figure 4.4: Impulse response of FIR filter [11] 
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4.2.2 simulation steps for m < 1 

Rayleigh distribution has been shown by many investigators, to be appropriate for 
modeling the fading statistics of a channel. However, it is possible that fading becomes 
more severe than Rayleigh fading as a result of high variability in the channel. The 
channel measurement results observed by Nakagami [11] indicate that the Nakagami- 
m distribution with m in the range of 0.5 < m < 1 is useful for modeling the fading 
characteristic of a channel when the fading is more severe than Rayleigh fading. 

Let z(t) = be a wide-sense stationary (WSS) complex random process that 

is characterized by the autocorrelation function 

R, (At) = E[z(t)z^ (t + At)] (4.4) 

and satisfies the following properties. 

1. r(t) is Nakagami-m distributed with the second moment f2 = ^^(O} and the 
fading parameter m limited in the range m G [0.5, 1). 

2. ^(t) is uniformly distributed over [0, 27r). 

3. r(t) and &(t) are mutually independent. 

Then z(t) can be represented by [15] 

z(t) = fj,(t)w(t) (4.5) 

where fi(t) is a WSS nonnegative random process and fj,^(t) follows a standard beta 
distribution with parameters m and 1-m, and w(t), being independent of /r(t), is a 
zero-mean WSS complex Gaussian process with an autocorrelation function 

R^(At) = E[w(t)w*(t + At)] (4.6) 

satisfying Ru,(0) = Q/m. Rx(At) and Rtp(At) are related by 

R,iAt) = R^(At).R^(At) (4.7) 

Where i2/!i(0) = m. The representation of z{t) by 4.5 is the desired simulation model. 
The simulation steps for the above model are [15]: 
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• /i(i) can be generated by where, J/(^) is a zero-mean unit 

variance Gaussian process with autocorrelation function p{At) = Ely(t)y{t+At)], 

is the CDF of y(t}, and F(x) = Ixi(m, 1 — m), 0 <x < 1, where, I^{a, b) 
is the incomplete beta function. y(t) is generated by using Rice’s sum of sinusoids 
[16]. 

• w{t) is generated independently by separately generating in-phase and quadrature- 
phase components based on the sum of sinusoids method. 

• z(t) is then generated by multiplying /i(t) and w(t) as in Eqn. 4.5 



Figure 4.5: Simulation model for Nakagami-m fading channel 
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4.3 Receiver 

The Nakagami-m channel model has been used to generate the raxnved signal given 
by Eqn. (3.1) for multipath propagation ss per the channel shown on Fig. 4.5. In the 
figure a is the attenuation coefficient and n is the time delay for path. The number 
of multipaths are taken to be 3 (same as IS-95 standard). 'Ihe receiver block diagram 
is shown in Fig. 3.2. 

The steps followed for simulation are; 

• The received signal is carrier demodulated and then pas,sed through the matched 
filter. The transfer function of the matche<i filt<‘r i.s complex conjugate of the 
filter used in transmitter. Since the filter response is symmetrical and real the 
response of matched filter is also same. Hence the same coefficients as shown in 
Table 4.1 are used in the receiver. 

• The matched filter output is then fed to a denuHiulator which is combination of 
L parallel demodulators ( L=3 in this casf;). 

• It has been assumed that the multipath delay.s art! known to the receiver, hence 
each demodulator despreads the received signals of each miillipath independently 
by correlating it with the 1 and Q PN .se<pience.s. Then the desj>read signal is 
averaged over a varying number of chips (/V,,) to e.stimate the relative path values 
(i.e. aicos<j>i and aisin^i) on the I-plmse and Q-phiise branche.H of Fig. 3.2. The 
value of Np should be Ies.s than or etpial to the mimlx'r of chips over which the 
attenuation and phase remains constant. 

• The estimated path values are inultiplied with lh«‘ ttrapread signal. Then the 
correlation with the Walsh seqmuice is carric'd out. 'Phe resulting signal from 
I-phase and Q-phase are added followed by summing the succeasive chips over a 
symbol duration of 04 chips . 
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In this simulation the performance of the forward link is studied over varying channel. 
This change in channel has been affected by changing the value of fading parameter m. 



Chapter 5 


Results and Conclusion 


Simulation of forward link of IS-95 standard has been carried out to study its perfor- 
mance over Nakagami-m channel. The number of resolvable multipaths is taken to be 3 
as specified in IS-95 standard. Nakagami-m channel model is considered to be a better 
model for mobile channel conditions, as it includes Rayleigh fading channel (m=l), 
one-sided Gaussian channel (m=0.5) and in the limit m ->■ +oo, non fading AWGN 
channels as special cases. Therefore, it is expected that the performance curves ob- 
tained will be more general. The transmitter and receiver have been implemented as 
per the block diagrams given in Fig. 3.1 and Fig. 3.2 respectively. Plots have been 
obtained for bit error rate (BER) vs. signal to noise ratio (SNR) under various channel 
conditions. 

In this work zero padded PN codes have been used as was discussed in section 4.1. 
Fig. 5.1 shows that the autocorrelation function is not much effected by padding a 
zero in the PN sequence. Also it is shown in Fig. 5.2 that zero padding does effects 
the autocorrelation function if the total length of PN sequence is less. 

The observations from results obtained by the simulations are as follow: 

• Case 1: m < 1 

A channel with more obstructions, such as in urban areas, can be more appro- 
priately modeled with Nakagami-m channel with m<l as compared to Rayleigh 
channel, which has been generally used to model such channels. Fig. 5.3 shows 
the BER vs. SNR curves that have been plotted for m = 0.5, 0.7 and 1.0 respec- 
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Figure 5.2: Autocorrelation function PN code of length 15. 




BER 
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tively. It also shows the Chernoff bound (Eqn. 3.14) for BER for m = 0.5. It can 



Figure 5.3: Simulation results for m= 0.5, 0.7 and 1.0 

be clearly seen that the curves are well below the bound and the increase in value 
of m results in decrease in BER. This agrees with the notion of Nakagami-m dis- 
tribution that the increase in value of m reduces the variance. It can also be seen 
that for m<l the performance is not as good as m=l (Rayleigh fading channel). 
This is again in conformity with the fact that Nakagami-m channel with m<l is 
worse when compared to Rayleigh fading channel. 

Case 2: m > 1 

Nakagami-m model with m >1 is a better model for suburban/rural areas, as 
the increase in value of m implies that the number of scatterers in the medium 
is reduced. Fig. 5.4 shows curves for different values of m = 1, 3, 5 and 9. The 
curves agree with the fact that with increase in m the performance gets better. 

• In Fig. 5.5 the curves for varying number of users, for the case when m = 3, 
have been plotted. It can be observed that as the number of users increase BER 
alf?n inrrfiases. But this increase is nominal because of the orthogonality between 








Chapter 1 


Introduction 


With the spectacular growth in the coininnuication facilities, the ileinaiKl <>1 liigli hat.i iat('| 
is increasing day Iw day. Most of tin' coniniunicaf ion f (S'linitjties liiice used sittgli' ciuricri 
modulation to transmit the data, which nec'ds coniplicaf ed e(juali/ei liluck.s at the n'ci'ivi'i’! 
to nullify the impairment made by channel. Though almost no comnninii'al ion .system rtiiis, 
without equalizers, hut we always try to minimize' tlu'ir m'e'd as the\ ate higlilv coiuidicatedi 
blocks. Normally more is the data raten uion^ comi>!icat<’d the ('i]itali/ei is. ilence cuiik'hI 
the idea of transmitting high speed data, (xtcujtying larg(' bandwidt h b\' dividing into sev(‘nil 
low speed data strea,ms occupying smaller bandwidt, In To do .so, vvi- tntnsniil iniinv ilatai 

I 

1 

streams on frequency division nmltiitlexed stibca,rriers. Ht'tice tin' equali/er nei'ded becoiiit's 
less complicated and sometimes it is comph'tely ('liininatt'd. Tin' s\ stems ba.sed on this idea 
are put under the fa,niil\' of multica.rri(‘r .syst('nis. 

1.1 Few implementations of multicarrier systems 

The principle of multicarrier was originally applittd in (’ollins' Kineplex Systems in ( !|. Then’ 
are claims of even early implementations. It has since' bee'ii called 'b\' many names and tisi'd 
with varying degree of sticcess in different nu’dia,. e.g.[cS] 

• FDM telephony group band modem. 

• Telephony voicee band modem. 

• Upstream cable mode’in. 


1 



users specific codes. 


• Fig. 5.6 shows the ChernofF bound and simulation result for Rayleigh fading 
channel. This plot shows that as expected the BER is always less than this 
upper bound. 

• The curves corresponding to upper bound given in Eqn. 3.17, ChernofF bound 
(Eqn. 3.14) and BER for m = 3 are shown in Fig 5.7. It can be seen that 
the upper bound is looser than ChernofF bound. This could be the result of 
the fact that the upper bound on BER given by Eqn. 3.16 is derived from the 
maximum value of the intgrand obtained over the interval of integration [13], 
whereas, ChernofF bound is derived from the density function. 


5.1 Scope for Future work 

This work has been carried out with the assumption that the number of paths and 
their delays do not change during call duration. It has also been assumed that delays 
are known to the receiver. In practical situations, number of paths and their delays 
keep on changing in accordance with the change in channel conditions. The work can 
be extended for practical situations by incorporating a timing acquisition circuit which 
adaptively estimates the number of paths and delays. Also the shape factor has been 
taken to be same for all the paths. Different values of rn for different ])aths may provide 
a more realistic performance. 
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